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1 = 2 pion-pion scattering phase shift in the continuum hmit calculated 
with two-flavor full QCD * 
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We present a calculation of the scattering phase shift for the 1 = 2 S-wave pion-pion system in the continuum 
limit with two-flavor full QCD. Calculations are made at three lattice spacings, using the finite volume method 
of Liischer in the center of mass frame, and its extension to the laboratory frame. 



1. Introduction 

Pilot studies of the scattering phase shift were 
carried out for the 1=2 S-wave two-pion system 
in Ref . |1I2| . These studies were made only at one 
lattice spacing employing the quenched approxi- 
mation. In this report, we present a calculation 
including two flavors of dynamical quarks, so that 
we can get rid of errors from unitarity violation, 
at three different lattice spacings to carry out ex- 
trapolation to the continuum limit. 

Our study is based on the finite volume method 
of Liischer formulated for the center of mass (CM) 
frame jj], and its extension to the laboratory 
frame by Rummukainen and Gottlieb While 
the inclusion of the lab frame adds little extra 
computational cost, the use of the lab frame al- 
lows significantly more dense sampling of the en- 
ergy states. 

We use the full QCD configurations [H| pre- 
viously generated with an RG-improved gauge 
action and a mean-field improved clover quark 
action at three gauge couplings /3 = 1.80, 1.95 
and 2.10, corresponding to a^^ « 0.92,1.3 and 
1.8 GeV, on 12^ x 24, 16^ x 32 and 24^ x 48 



lattices. The physical volumes are about 2.5 
fm. Four hopping parameters corresponding to 
m^r ~ 1.0,0.9,0.75 and 0.55 GeV arc taken at 
each (3 for chiral extrapolation. The numbers of 
the configurations vary from 380 to 720. 

2. Methods 

On a finite volume of L'^ the interaction be- 
tween two pions shifts the two-pion energy E. 
In the CM frame, writing E = 2^y^r^ + J? 
with = (27r/L)^n (n being non-integer), 
the phase shift 6{p) is given by the relation 
tan(^(p)) = TT^^^^/^f Zooi^',n) where Zooii',n) = 
{1/^/4tt) X^K'^ ~ ^)~% ' being integer vectors. 

In the laboratory frame, the two-pion energy 
El and the sum of the two pion momenta P = 
Pi + P2 are given by = (27r/L)^TO = — 
p2)/4- 
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- m^. The phase shift is now tan((5(p)) = 
77r"/'V^/^(^(l; m) where 7 = ElI \/E\~ 

and ZP^{i;m) = {l/V^)J2rir^-^)''- Here 
the summation is taken over the vectors r that 
are constructed as r = J^^H + d/2) where d = 
LP/27r. 

The pion four-point function behaves as a sum 
of exponentials due to the presence of a num- 
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Figure 1. Ratio of the four-point to the square of 
two-point functions. The values before diaginal- 
ization Rn{t) are compared to those after diagi- 
nalization Wn(t). Note that Wo{t) is completely 
overlayed on Ro{t). 

ber of states having the same quantum num- 
bers. To resolve the energy eigenvalues, we con- 
struct in each frame the pion four-point func- 
tion matrix Cnm{t) = {0\fln{t)flm{ts)\0) , where 
ts is the source point and Jl„(i) = 7r(pi)7r(p2) 
is the two-pion operator with momenta pi that 
depends on the state n. The energy eigenvalues 
An(t) = exp[—En{t — to)] are obtained by diago- 
nalizing the matrix C-^/'^{to)C{t)C-^/'^{to) with 
to a reference time [HI . A cut-off N is introduced 
in the number of the energy states considered. 

We work with the CM and two laboratory 
frames, LI and L2. The sums of the two pion mo- 
menta P in LI and L2 are Pli = (1, 0, 0) x (27r/L) 
and Pl2 = (1,1,0) X {2Tr/L). We choose = 3 
in the CM frame and 4 in the lab frame. In each 
frame, 3 x 3 or 4 x 4 matrices are diagonalized and 
the phase shift is determined for the energies of 
the ground (n = 0) and the first excited {n — 1) 
states. 

FigureHillustrates the effect of diagonalization. 
We compare the ratio Rn{t) = C„„(i)/G^j(t), 
where G„(t) = {0\n{pn,t)n{-pn,ts)\0) withp^ = 
(27r/L)^n, before the diagonalization (denoted by 
open symbols) to Wnit) = \n{t)-Gl{to)/Gl{t) af- 
ter the diagonalization (closed symbols) normal- 
ized a,i t = ts- For the ground state, Wo{t) ~ 
Ro(t) so that the diagonalization is unnecessary. 
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Figure 2. Scattering amplitudes defined by (1) 
for the heviest and lightest pion masses at (3 = 
1.80 (top panels) and 2.10 (bottom panels). The 
dashed lines show the global polynomial fit. 

For the n = 1 state, however, the ratio shows 
an exponential fall-off only after diagonalization: 
Ri{t) blows up at large t. We apply an exponen- 
tial fit to Wi{t). The energy eigenvalues for the 
n — and 1 states change very little even if we 
adopt the cut-off N = 2 instead of 3 and 4. 

3. Phase shifts 

We define the 'scattering amplitude' 

Tip, m^, a) ^ ta,n{d {p))E^/p (1) 

where Et^ = \/m^ + p^. This gives the scat- 
tering length flo in the zero-momentum limit: 
limp^o T{v, m^, a) = aom^. 
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Figure 3. Phase shifts at the physical pion mass 
for the three values of j3. 

Figure El shows the scattering amplitude at 
/3 = 1.80 and 2.10 for the lightest and heaviest 
pion masses. The CM„ refers to the amplitude 
obtained from the n-th state in the CM frame, 
and the Ll„ and L2„ to those in the lab frames. 
Note that we have three points between the two 
CM data, showing the improvement of the mo- 
mentum sampling. 

We carry out a global fit to the T at each (3 with 
a polynomial T(p, TO^,a) — Aiom^ + A2om^ + 
Asaml + Aqip^ + Aiimlp^ + A2imlp'^ . One data 
point (shown with open symbol) is excluded from 
the fit, since the Wn{t) does not show a good 
exponential behavior. We find a reasonable fit to 
the data for both CM and lab frames. 

In Figl^l the phase shift at the physical pion 
mass = 0.14 GeV, as calculated from the fit, 
decreases with /3, showing the presence of an 0(a) 
effect. 

The continuum limit is taken by extrapolating 
the fits linear in a to the scattering amplitude 
at the physical pion mass for each momentum, 
T(j}, TO^, a) = T{p, TO^) + aT°^{p, m^). The phase 
shift in the continuum limit is presented in Fig01 
with the dashed line, associated by a band of error 
bars. This is compared to the experimental data. 

The scattering length obtained by taking the 
zero-momentum limit of T{p, mjr) in the contin- 
uum, aorriTr = —0.0488(49), may be compared 
with the prediction of chiral perturbation theory: 
aom^ = -0.0444(10) |7|. 




Figure 4. Phase shifts in the continuum limit (in- 
dicated by the dashed curve) . The band of bars 
shows our estimate of errors. The solid curve is 
parametrized by experimental inputs |2| and the 
symbols are the experimental data. 



This work is supported in part by Grants-in- 
Aid of the Ministry of Education (Nos. 12304011, 
12640253, 13135204, 13640259, 13640260, 
14046202, 14740173, 15204015, 15540251, 
15540279, 15740134 ). The numerical calcu- 
lations have been carried out on the parallel 
computer CP-PACS. 

REFERENCES 

1. H. R. Fiebig et ai, Few-Body Syst. 29 (2000) 
95. 

2. CP-PACS Collaboration: S. Aoki et at, Phys. 
Rev. D67 (2003) 014502. 

3. CP-PACS Collaboration: A. Ah Khan et ai, 
Phys. Rev. Lett. 85 (2000) 4674; Phys. Rev. 
D65 (2002) 054505. 

4. M. Liischer, Commun. Math. Phys. 105 
(1986) 153. 

5. K. Rummukainen and S. Gottlieb, Nucl. 
Phys. B450 (1995) 397. 

6. M. Liischer and U. Wolff, Nucl. Phys. B339 
(1990) 222. 

7. G. Colangelo et ai, Nucl. Phys. B603 (2001) 
125. 



